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^ ■ ABSTRACT 

I ^ I ' Aims. We report the detection of transit timing variations (TTVs) well in excess of one hour in 

the Kepler multi-planet candidate system KOI 806. This system exhibits transits consistent with 
^ ■ three separate planets - a Super-Earth, a Jupiter, and a Saturn - lying very nearly in a 1:2:5 reso- 

, nance, respectively. 

Methods. We used the Kepler public data archive and observations with the Gran Telescopio de 
Canarias to compile the necessary photometry. 

Results. For the largest candidate planet (KOI 806.02) in this system, we detected a large transit 
timing variation of -103.5±6.9 minutes against previously published ephemeris. We did not ob- 
tain a strong detection of a transit color signature consistent with a planet-sized object; however, 
we did not detect a color difference in transit depth, either. 
■ Conclusions. The large TTV is consistent with theoretical predictions that exoplanets in res- 

onance can produce large transit timing variations, particularly if the orbits are eccentric. The 
presence of large TTVs among the bodies in this systems indicates that KO1806 is very likely to 
be a planetary system. This is supported by the lack of a strong color dependence in the transit 
depth, which would suggest a blended eclipsing binary. 

Key words. Techniques: photometric ~ planetary systems - star: individual: KOI 806, KIC 
3832474 

1. Introduction 



Transit timing variations (TTVs) are proving to be a very valuable tool in exoplanet research. In 
systems with transits of different periodicities, they may serve to constrain the mass of the transiting 
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Fig. 1. GTC photometry of a transit of candidate KOI 806.02 in g' (top) and z' (bottom). The solid 
line is the fitted transit center and the dashed line its 3cr confidence limit, the shaded region defines 
the transit center +3cr predicted by Bll, and the dotted hne is mid-way through the egress. The 
time of transit center predicted by the ephemeris of Bl 1 is 2455740.5146 + 4.2 x lO""', about l.lh 
earlier. The lower figure shows the observed g' - z' during the transit egress, along with a modeled 
color signature using limb darkening models from Claret & Bloemen (1201 11 1. The points shown are 
the mean of each 4 x 5i exposure sequence. 

bodies and thereby verify them as planets, without further observations by e.g. radial velocities. 
Furthermore, TTVs may lead to the discovery of further non-transiting planets. Indeed, TTVs have 
been instrumental in characterizing multiple systems found by Kepler, including Kepler-9 (Holman 
et al. l20T0l l. Kepler- 1 1 (Lissauer et al. l20T0l l. Kepler- 16 (Doyle et al. lSOTlT l. Kepler- 18 (Cochran et 
al. l201 11 1, and Kepler- 19 (Ballard et al. l201 11 1, where an additional, non-transiting planet has been 
identified solely through TTVs. 
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Fig. 2. ALFOSC/NOT spectram for KOI 806 (top) along with a K5 standard star observed with the 
same instrumental configuration on the same night. The two spectra are smoothed by 5 pixels and 
offset by 0.3 for clarity, with the pertinent atomic Unes indicated. 
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Fig. 3. Pubhc Kepler photometry for the two released transits of KOI 806.01 (top), KOI 806.02 
(bottom) along with the best fit transits. 
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Fig. 4. Transit timing variations for the KOI 806 candidates vs. the periods and Tqs given in Bll. 
All of the candidates in the system exhibit high significance (cr > 6) TTVs, with KOI 806.01 and 
KOI 806.03 doing so in just the Kepler Q0-Q3 data. KOI 806.01 has a min/max difference of over 
1.5 hours for just the first two transits, KOI 806.02 1.8 hours, and KOI 806.03 1.4 hours. These 
values are all tabulated in Table [T] 

Since its public release in Feb. 2011 of the list of Kepler planet candidates (Borucki et al. 
1201 II hereafter Bll), the system KOI 806 has been an obvious candidate for large TTVs. The 
innermost two planet candidates, designated 806.03 and 806.02, (hereafter '05' and '02') are given 
with periods of 29.1654+0.0012days and 60.32875±0.00037, respectively, and a further candidate, 
806.01 (hereafter '07'), is listed with P=143. 1814+0.0027 days. The periods of 02 and 03 lie very 
close to the 2:1 resonance, which can strongly enhance TTVs relative to a non-resonant state (e.g. 
Agol et al. I2005I I. However, the analysis of transit timing observations from Kepler observing 
segments QO - Q2 (corresponding to data up to JD 2455091.48) by Ford et al. (1201 11 1 stops short 
of listing this system among clear candidates for detected TTVs. The apparent reason is that in 
data from Q0-Q2, only the shortest-periodic 03 exhibited possible TTVs, whereas 01 and 02 had 
only one transit in these segments. Even so, the TTVs predicted by Ford et al. for all candidates in 
KOI 806 are among the highest of the candidates in Bll; Ford et al. predict that the min-to-max 
TTVs over 3.5 years could range from ~ 15 minutes for 07 to almost 5 hours for 03, the smallest 
candidate with the shortest orbit. 

2. Observations 

We used the 10.4m Gran Telescopio Canarias (GTC) to observe a transit of candidate 02 with 
semi-simultaneous time series photometry in g' and z', 2 widely separated colors chromatically. 
These observations were undertaken in an attempt to verify that the transiting body is indeed of 
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planetary size, as it should then show characteristic color variations due to the interplay between 
the relative size of the transiting body and differential stellar limb darkening (Tinglev f2004l l. While 
such detection of the color "signature" of an exoplanet can be difficult, 02 has a deep transit (almost 
2%) and a long transit duration (6.6 hrs), which made it an excellent candidate for the attempt. We 
obtained approximately 1800 images with GTC/OSIRIS in five hours on the night of the 27/28th of 
June, 2011, alternating four 5-second exposures in g' (468.8 nm) with four 5-second exposures in 
z' (893.1 nm). Subsequent analysis using optimized aperture photometry with the ' Vaphot' package 
(Deeg & Dovle l200TT l for IRAF showed that the transit unexpectedly ended 100 min earlier than 
anticipated by the period and duration given in Bll, a detection of a TTV with a significance of 
~ 15cr. Given the presence of 3 candidate transits in this system and the high significance of this 
TTV detection, it is extremely likely, just based on this, that KOI 806 harbors a planetary system 
of at least three planets. 

In the light of this very strong TTV detection, the original science goal based on multi-color 
photometry becomes relegated to supporting evidence. Despite the TTV, we captured the crucial 
part of the egress where the planet crosses the limb of the star, which results in a weak detection of 
2i g' - z' color signature (on the order of one millimag), consistent with the modeled one, though 
offset (Fig.[T]l. This offset is plausible, given the clear evidence of transparency variations just after 
the transit, revealed by increased scatter in the photometry. This could in turn affect the measured 
out-of-transit flux. 

However, even without a true detection of the transit color signature, we can constrain the the 
nature of any false positive. The photometry does clearly excludes the possibility of any strong de- 
pendency of transit depth on color, which one would expect from a blend with any color difference 
between the components, as proposed by Tingley (120041) and utilized by O' Donovan et al. ( I2006I I. 
O'Donovan et al. ( 120071 1. Cochran et al. (1201 11 1, and Ballard et al. (120111 ). We derived an equation 
that is similar in principle to BLENDER (Torres et al 120041 1201 11 1, but much simplified, to limit the 
possible color difference between the target star and any contaminating eclipsing binaries (CEBs). 
It relates the blended echpse depth in z' (dcEB,g) to that in g' using only the depth of the unblended 
ecUpse depth in g' and z' ( dEB,g and c/^b,;, respectively) and the g' - z! color difference between 
the target star and the blending eclipsing binary (A(g' - z'))- 

1 + dEB,z.fEB,z 

dcEB,z = \ 7 . (1) 

i + Teb,z 



where 

jEB,z - Jtarger,g^^ 



1 - dcEB,g 
dcEB.g - dEB4 



where we set ftarget.g - 1, which yields a differential blended eclipse depth in z' ■ If we set dEB,z - 
dEB,g, Mg' - z') must be less than about 0.02. Relaxing this constraint to allow dEB.z to vary as 
much as 10% from dEB,g, a reasonable value, constrains A(g' - z') to about 0.25 for lightly blended 
eclipses {fiargei ~ /eb) to 0.05 for highly blended eclipses {ftarget ^ Ieb)- While this does not 
absolutely rule out the possibility that the transit of 02 may be due to a blend, it does suggest that 
it is unlikely. 

While the GTC showed a clear advance of the transit egress, we would not have been able to 
detect any changes in transit duration due to the partial coverage of the transit. To verify the stellar 
parameters given in the Kepler Input Catalogue for the host star, KIC 3832474 (T^s - 52Q6K, 
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Table 1. Best-fit transit times for the planet candidates KOI 806.01, KOI 806.02 and KOI 806.03, 
and their offset (O-C) against the ephemeris given in B 11 . is the transit number and the dates are 
in BJD-2454900. 
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258.6832 ± 1.2 x 10^ 


-3 


-1.7 ± 10.6 



'From egress observed by GTC 



log^ = 4.53), we took spectra of KOI 806 and several other similar stars with the 2.5 meter Nordic 
Optical Telescope (NOT), using grism #5 at a resolving power R=830 over the wavelength interval 
5100-9000A. The spectra were calibrated in wavelength against Ne lines from arc images taken 
with the same instrumental configuration as the target and flux normalized to unity at the wave- 
length interval 7400-7500A. Uncertainty in instrumental response correction is less than 10%. We 
corrected the spectra for telluric absorption, though some residuals at the strongest telluric oxygen 
band at around 7600A remain due to different airmasses of target and telluric standard star (see 
Fig.|2]i. We obtained the spectral classification of the host star via atomic line identification (Mgl, 
Nal, Ha, and Call), which indicated a spectral type of G9-K0+2 and a slightly higher temperature 
(Teff ~ 5250K) than given in Bll. We did not have the necessary spectral precision to estimate 

This leaves the question of the radius of the star open, given by B 1 1 as = 0.88.Ro- A Jupiter- 
sized planet transiting a star of this size would cause a transit with a depth of about 1.4%, whereas 
the observed transits of 02 have a depth of 2.0%. Therefore, for the transits of 02 to be due to a 
planet, 02 needs to be ~ 20% larger than Jupiter or the star has to be smaller than expected by a 
similar amount - or, most likely, some combination of the two. 

3. Analysis of transits 

With the recent release of Kepler Q3 data, covering up to 16 Dec 2009 (JD 2455182.00), we were 
able to extend our analysis, as second transits of both 01 and 02 are now available. We performed 
fits (see discussion below) to the transits in Kepler data of each the three candidates (see Fig.[3]for 
01 and 02, 03 similar but not shown). 

Our parameter estimation process is based on a Bayesian methodology. We first search for the 
rough global maximum likelihood solution using the differential evolution (DE) global optimiza- 
tion method (Price, Storn & Lampinen 120051 Storn & Price lI99"7l l. and then derive the parameter 
posterior distributions using Markov Chain Monte Carlo (MCMC) (Ford l20()5l l. We parameterize 
the transits by the transit center Tc, squared radius ratio k, period p, impact parameter b, and re- 
ciprocal of the half duration of the transit T, similar to Kipping (120 IOD and Bakos et al. (120 1 01) . 
Further, we use the linear limb darkening law with a free limb darkening coefficient u. 
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We use unconstraining (uninformative) priors for the parameters estimated for the two largest 
candidates, 01 and 02. For the smallest candidate, 03, we use normal distribution priors to constrain 
the values allowed for b, Y and u. The prior centers for b and T are adopted from the released Kepler 
values and widths are chosen conservatively, as we prefer to overestimate the errors slightly rather 
than underestimate them. The limb darkening prior is derived from the limb darkening posteriors 
of the fits to 01 and 02. 

For the fits to the long-cadence Kepler data, we subsample our model by 10 (i.e., each final 
model light curve point is an average of 10 model points uniformly spread over the time of a single 
Kepler exposure). This is necessary since long exposures modify the transit shape (Kipping |2010l) . 
The Kepler observations and the corresponding fits can be seen in Fig. [3] for 01 and 02. 

The fit to our GTC observations of candidate 02 is done simultaneously in both colors. The only 
major difference to the single color Kepler data is that we have separate limb darkening coefficients 
for each color. Since we have only a partial transit, we cannot infer the transit duration. Instead, 
we use an informative prior of 5.7 ±0.1/; from the transit duration posteriors of the Kepler transit 
fits, and similarly, a prior for the impact parameter ofb- 0.550 + 0.075. The rest of the parameters 
have uninformative priors for conservative error estimates. The resulting residuals have a cr^ms of 
1 .9 millimags in g' and 2. 1 millimags in z' per 4x5 exposure sequence. 

With the time of the transit center being one of the free fit parameters, we were able to measure 
the periods from the first two transits of 01 and 02 and the first seven transits of 03, respectively, for 
comparison to the period given by Bl 1. In combination with our own observations with the GTC, 
we can confirm that large and statistically very significant TTVs are apparent in this system (Fig.|4] 
and Table 1). 

4. Conclusions 

We have detected strong TTVs in the KOI 806 system, affecting the timing of the transit far more 
than what we had expected based on the uncertainties indicated by the ephemeris of Bll. The 
magnitude of the TTVs seen in this system are in fact not surprising, as Holman & Murray (I2005I I 
have argued that TTVs may scale in a roughly linear fashion with period, though TTVs are highly 
degenerate and strongly dependent on initial conditions (Nesvorny & Morbidelli l2008l Veras, Ford 
& Pavne i20 irn . If this linear relationship holds, then the observed TTVs in this system are similar in 
scale to those observed in, for example, the Kepler- 11 system (Lissauer et al. 120 101) - with TTVs of 
tens of minutes at periods of 10 - 47 t/. Our observations of 02 indicated a TTV of - 103.3 + 6.9 min 
relative to the transit ephemeris given by the Kepler candidate list (Bll). These values are based on 
a linear fit to all the transit that occurred during Kepler quarters Q0-Q5. Over this period of time, 
02 transited 5 times, yet the Bll values imply a transit ephemeris error of only ~ 4 min for the 12th 
transit, the one we observed with GTC (see FiglU. We performed an independent analysis of the 
publicly available transits (Q0-Q2), finding significant differences to the periods and Tos given in 
Bll. For example, we derive a period of 60.33465 + 0.00042 d for 02, which results in a TTV for 
the 12th transit of ~ 3.2 /i. However, as much of the Q0-Q5 data are still proprietary, we cannot 
perform a full, independent revision of the method and results that led to the ephemeris given by 
B 1 1 . A TTV well in excess of one hour is certain, though, regardless of which period is used, 
either ours from the first two transits or that from Bll. The second result of our GTC observations 
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is the absence of any color dependency of the transit depths, the presence of which would reject 
the candidate, along with a color signature during egress that suggests that an object with a radii 
consistent with a planet is causing the transit, albeit with a low statistical significance. 

While our data do not allow for a detailed analysis of the system, we consider that the exis- 
tence of a very strong TTV and the absence of a color dependency of the transit depth strongly 
supports that the transits in KOI 806 are caused by planetary companions. Two possible scenarios 
for massive companions exist: either the planet is orbiting one star or is circumbinary. Transiting 
circumbinary planets have a complex pattern of transits due to the relatively large motion of the 
stellar components as they orbit around the center of mass, which is clearly not seen is this system 
(Deeg et al. 12000 1 Doyle et al. 1201 11 1. TTVs in systems where the planet orbits one star, called 
S-type orbits, are caused by the light time effect (e.g. Kopal [19591 as distance between the transited 
star and the observer changes as the stars orbit around the center of mass. While this can cause 
large TTVs over time, a 1.7 h TTV requires an acceleration amounting to a change of about 13 
AUs. Such a change could only be explained by some exotic, invisible and very massive object 
(e.g. a black hole), as a small, faint star could not move the primary star so far in the less than 3 
years over which the TTVs have been measured, and is therefore very unlikely. Moreover, it could 
still not explain the TTVs measured in the other planets in the system, assuming the values in B 1 1 
are correct. 

Given the long periods of the putative planets in KOI 806, it is entirely possible that this system 
will not be fully characterized by the end of the Kepler mission - particularly if evidence of other 
unknown, perhaps longer period planets arises. However, the deep transits and large TTVs make it 
possible to monitor this system from ground over the long term - an endeavor likely to be scien- 
tifically worthwhile, depending on the details that emerge from the analysis of the data currently 
being acquired by the Kepler mission. 

A final aspect of unknown TTVs in this and other planetary systems pertains to photometric 
follow-up. It appears that the transit ephemeris given in B 1 1 do not (or cannot) necessarily account 
for TTVs in all cases, resulting in ephemeris errors that significantly underestimate the true timing 
error Therefore, observers intending to follow up Kepler candidates, particularly in known multi- 
candidate systems where large TTVs may be present, need to include this fact when planning their 
observations. 
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